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Abstract. In this paper, firstly, we study analytically the topological features of a 
family of hierarchical lattices (HLs) from the view point of complex networks. We 
derive some basic properties of HLs controlled by a parameter q: scale-free degree 
distribution with exponent 7 = 2 + null clustering coefficient, power-law behavior 
of grid coefficient, exponential growth of average path length (non-small- world), fractal 
scaling with dimension ds = ^wf^ , and disassortativity. Our results show that 
scale-free networks are not always small-world, and support the conjecture that self- 
similar scale-free networks are not assortative. Secondly, we define a deterministic 
family of graphs called small- world hierarchical lattices (SWHLs). Our construction 
preserves the structure of hierarchical lattices, including its degree distribution, fractal 
architecture, clustering coefficient, while the small-world phenomenon arises. Finally, 
the dynamical processes of intentional attacks and collective synchronization arc 
studied and the comparisons between HLs and Barabasi- Albert (BA) networks as well 
as SWHLs are shown. We find that the self-similar property of HLs and SWHLs 
significantly increases the robustness of such networks against targeted damage on 
hubs, as compared to the very vulnerable non fractal BA networks, and that HLs have 
poorer synchronizability than their counterparts SWHLs and BA networks. We show 
that degree distribution of scale-free networks does not suffice to characterize their 
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synchronizability, and that networks with smaller average path length are not always 
easier to synchronize. 

PACS numbers: 89. 75. Da Systems obeying scaling laws -05.45.Df Fractals -36.40.Qv 
Stability and fragmentation of clusters -05.45.Xt Synchronization; coupled oscillators 
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1. Introduction 

Topological characteristics, such as scale-free degree distribution, small-world effect, 
fractal scaling and degree correlations, have recently attracted much attention in 
network science. The last few years have witnessed a tremendous activity devoted 
to the characterization and understanding of networked systems [U EJ [3j Hj. Small- 
world property [5] and scale-free behavior [6j are two unifying concepts constituting 
our basic understanding of the organization of real-life complex systems. Small-world 
property refers to the one that the expected number of edges (links) needed to pass 
from one arbitrarily selected node (vertex) to another one is low, which grows at most 
logarithmically with the number of nodes. Scale-free behavior means the majority of 
nodes in a network have only a few connections to other nodes, whereas some nodes are 
connected to many other nodes in the network. This poses a fundamental question how 
these two characteristics are related. It has been observed that small-world property 
and scale-free behavior are not independent p]: scale-free networks, normally, have 
extremely short average path length (APL), scaling logarithmically or slower with 
system size. Is this universal? 

In fact, the above mentioned two properties (i.e. small-world property and scale- 
free behavior) do not provide sufficient characterizations of the real-world systems. It 
has been observed that real networks exhibit ubiquitous degree correlations among 
their nodes [H [9j [TOj EH 121 E31 ITS] . This translates in the observation that the 
degrees of nearest neighbor nodes are not statistically independent but mutually 
correlated in practically every network imaginable. Correlations play an important 
role in the characterization of network topology, and have led to a first classification 
of complex networks A series of recent measurements indicate social networks 

are all assortative, while all biological and technological networks disassortative: in 
social networks there is a tendency for the hubs to be linked together, in biological and 
technological network the hubs show the opposite tendency, being primarily connected 
to less connected nodes. Correlations are now a very relevant issue, especially in view 
of the important consequences that they can have on dynamical processes taking place 
on networks [151 ESI EH EE] • 

Recently, it has been discovered by the application of a renormalization procedure 
that diverse real networks, such as the WWW, protein-protein interaction networks 
and metabolic networks, exhibit fractal scaling and topological self-similarity [T9J, EH1 
I2TI [22| [23] . Fractal scaling implies that in a network the minimum number of node- 
covering boxes Nb of linear size £b scales with respect to ^ as a power-law Nb, with an 
exponent that is given by a finite fractal dimension ds [19] . Self-similarity refers to the 
invariant scale-free distribution probability to find a node with degree k, P(k) ~ fc -7 , 
i.e. the exponent 7 remains the same under the renormalization with different box 
sizes [T9l |21] . In complex networks, fractality and self-similarity do not always imply 
each other: a fractal network model is self-similar, while a self-similar network is not 
always fractal [23]. One can obtain the fractal dimension ds by measuring the ratio of 
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Nb over the total number of nodes N in the network, which satisfies Nb/N ~ Iff 11 1 . After 
renormalizing the networks, the degree /cb(^b) of a node in the renormalized network 
versus the largest degree khub inside the box that was contracted to one node with degree 
khub in the renormalization process exhibits a scaling behavior: kg (&b) = s(£B)kh u b, 
where s(£b) is assumed to scales like s(£b) ~ £~B dk with dk being the degree exponent of 
the boxes. In self-similar scale-free networks, the relation between the three indexes 7, 
ds and dk satisfies 7 = 1 + d B /d k [19|. 

Correlations and topological fractality are important properties for many real- 
world complex systems. Then a natural fundamental question is raised how the 
two characteristics relate to each other. Recent researches [201 122] have shown that 
self-similar scale-free networks are not assortative, and the qualitative feature of 
disassortativity is scale-invariant under renormalization. Moreover, self-similarity and 
disassortativity of scale-free networks make such networks more robust against a sinister 
attack on nodes with large degree, as compared to the very vulnerable non fractal scale- 
free networks [2D] . 

However, do small-world property and scale-free behavior always go along? How 
do systems have evolved into self-similar disassortative scale-free networks? How the 
dynamical processes such as intentional attack and synchronization are influenced by 
the topological fractality and disassortativity of scale-free networks? Such a series of 
important questions still remain open. To relate these questions, in this paper we 
therefore launch a study seeking a better understanding of the relations among these 
topological properties. 

It is of interest to study above important questions with deterministic methods. 
Because of their strong advantages, deterministic network models have received much 
attention [251 12H1 [2TI [28l [291 EQl EU [32l [331 EH [35l EHl EZl EH [39] . First, the method of 
generating deterministic networks makes it easier to gain a visual understanding of how 
networks are shaped, and how do different nodes relate to each other [25]; moreover, 
deterministic networks allow to compute analytically their topological properties, which 
have played a significant role, both in terms of explicit results and a guide to and a 
test of simulated and approximate methods I251I2S1I2ZII2E1I2S1ED1EI1E21E31E31ES1 
I3P1 [37] [38] [3D] . On the other hand, deterministic networks can be easily extended to 
produce random variants which exhibit the classical characteristics of many real-life 
systems 0D1 SB S2J M, IB] • 

Inspired by the above mentioned questions, here we first introduce a deterministic 
family of networks. These networks are called hierarchical lattices (HLs), which yield 
exact renormalization-group solutions [45j |46l HTJ EHJ |49l [50]. From the perspective 
of complex network, we show that HLs are simultaneously scale-free, self-similar and 
disassortative, but lack the small-world property. Then we present a deterministic 
construction of a class of small- world hierarchical lattices (SWHLs), which preserve 
the basic structure properties including power law degree distribution, self-similarity, 
and disassortativeness, while lead to the small- world effect. Finally, we investigate the 
effects of network structures on the dynamics taking place in them. 
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Figure 1. Iterative construction method of the hierarchical lattices for some limiting 
cases. 

This article is organized as follows. In Sec. II, we introduce the construction of 
hierarchical lattices (HLs) and study their topological features including the degree 
distribution, moments, clustering coefficient, grid coefficient, self-similarity, degree 
correlations, and average path length (APL). The detailed exact derivation about APL 
is shown in Appendix. In Sec. Ill, we propose the deterministic construction of the 
small-world hierarchical lattices (SWHLs) and study their properties. In Sec. IV, attack 
tolerance of HLs is studied and the comparisons between HLs and Barabasi- Albert (BA) 
networks are shown. In Sec. V, we do a comparative investigation of synchronization in 
HLs, SWHLs and BA networks. Sec. VI is devoted to our conclusions. 

2. Hierarchical lattices 

In this section, from topological perspective of complex networks, we present the 
construction and the basic properties such as degree distribution, clustering coefficient, 
average path length (APL), fractality, and correlations of the hierarchical lattices (HLs). 

2.1. Construction of the Lattice 

The hierarchical lattices [H] are constructed in an iterative manner as shown in Fig. [TJ 
We denote the hierarchical lattices (networks) after t generations by H(q,t), q>2 and 
t > 0. The networks are constructed as follows: for t = 0, H(q, 0) is an edge connecting 
two points. For t > 1, H(q,t) is obtained from H(q,t — 1). We replace each of the 
existing edges in H(q,t — 1) by the connected cluster of edges on the right of Fig. [U 
The growing process is repeated t times, with the infinite lattices obtained in the limit 
t — > oo. Figure [2] shows the growing process of the networks for three particular cases 
of q = 2, q = 3, and q — 4. It should be noted that in the hierarchical lattice of q = 2 
case [15] , the Migdal-Kadanoff J5TJ [52] recursion relations with dimension 2 and length 
rescaling factor 2 are exact. 

Griffiths and Kaufman provided two explanations for the construction of the 



Self- similarity, small-world, scale-free scaling, disassortativity, and robustness in hierarchical lattices^ 





a) 

q=2 




t=0 



t=1 



t=2 



Figure 2. Examples of the hierarchical lattices for some particular case of q = 2, 
q = 3, and q = 4, showing the first three steps of the iterative process. 



hierarchical lattices [37], which are called "aggregation" and "miniaturization". In 
essence, these two interpretations reflect the self-similar structure of the hierarchical 
lattices, which allow one to calculate analytically their topological characteristics. 

Next we compute the numbers of nodes (vertices) and links (edges) in H(q,t). Let 
L v {t) and L e (t) be the numbers of vertices and edges created at step t, respectively. 
Note that each of the existing edges yields q nodes, and the addition of each new node 
leads to two new edges. By construction, for t > 1, we have 

L v {t) = qL e {t - 1) (1) 

and 



L e (t) = 2L v (t). 

Considering the initial condition L v (0) = 2 and L e (0) = 1, it follows that 
L v (t)=q(2qy- 1 

and 

L e (t) = {2qf. 

Thus the number of total nodes N t and edges E t present at step t is 



(2) 
(3) 
(4) 

(5) 
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and 

Et = L e (t) = (2qf, (6) 

respectively. 

2.2. Degree distribution 

Let ki(t) be the degree of node i at step t. Then by construction, it is not difficult to 
find following relation: 

ki(t) = qki(t-l), (7) 

which expresses a preference attachment [Sj. If node i is added to the network at step 
ti, ki(ti) = 2 and hence 

hit) = 2^. (8) 

Therefore, the degree spectrum of the network is discrete. It follows that the degree 
distribution is given by 

L v (0) 2 



P(k) 



Nt - g(2g)«+3g-2 f ° r *i - 

L V (U)_ qjqh fQTt . >1 (9) 

2q-l 

otherwise 



and that the cumulative degree distribution [HI [26] is 

p (k) - V Lv{p) - q{2q)U + 39 " 2 (10) 

In * 

Substituting for ^ in this expression using ti — t— gives 

_ g(2g) i (|)"^+3g-2 
cuml j ~ g(2g)* + 3q - 2 

for large t. (11) 



, _f1 + ln2 



2 

So the degree distribution follows a behavior power law with the exponent 7 = 2 + 



ln2 
In g ' 

For q = 2, Eq. ( TTTT) recovers the result of the particular case p = previously obtained 
in Ref. [26] . 



£5. Moments 

Information on how the degree is distributed among the nodes of a undirected network 
can be obtained either by the degree distribution P(k), or by the calculation of the 
moments of the distribution. The n-moment of P(k) is defined as: 

(k n ) = J2k n P(k). (12) 
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The first moment (k) is the mean degree. At arbitrary step t, the average vertex 
degree of H(q, t) is 

2E t 2(2g-l)(2 g )' 

<*>' = 74 = g (2,)- + 3,-2 ' (13) 

For large t, it is small and approximately equal to a finite value 4 — -. 

We can also calculate higher moments of the distribution P(k). For instance, the 
second moment, which measures the fluctuations of the connectivity distribution, is 
given by 

(k 2 )t = ^Y,n v {U)[k{t u t)}\ (14) 

* u=o 

where k(ti, t) is the degree of a node at step t which was generated at step £j. This 
quality expresses the average of degree square over all nodes in the network. It has large 
impact on the behavior of dynamical processes taking place in networks 
Substituting Eqs. d2J), ((SD and (jSJ) into Eq. ffT4|) . we derive 

4g(2g - 1) q* - 2* 2(2g - l)g* 



^ 4 f (3t + 3) 
4* + 2 
4(g-l)(2g-l) ,q 



for q > 2 
for g = 2 

— . . i for g > 2 . , 

g(g-2) \2) H (15) 

3(t + l) forg = 2 

In this way, second moment of degree distribution (k 2 ) has been calculated explicitly, 
and the result shows that it becomes infinite for large t. In fact, because the degree 
exponent 7 < 3, all n-moments (n > 2) diverge. For example, we can analogously get 
the third moment as 

2*+2^3t+3 _ 2 t + 1 q 3t + 2 _|_ 2*+ 3 g 3 *+ 1 _ 2 t+2 q 3t — 2 2t+A q t+1 + 2 2t+3 q t 
\ k )* = 2 2 *g i+3 - 2 2t+1 q t+1 + 3 ■ 2*g 3 - 3g ■ 2* +1 - 2* +1 g 2 + 2 t + 2 

For the special case of g = 2, it reduces to 
q . 9 3t _ a . ? 2t 

= 22t+ \ , (17) 
which diverges as an exponential law when t is very large. 



2.4- Clustering Coefficient 

The clustering coefficient defines a measure of the level of cohesiveness around any given 
node. By definition, the clustering coefficient |5] Cj of node i is the ratio between the 
number of edges e« that actually exist among the ki neighbors of node i (i.e. the number 
of triangles attached to a vertex %) and its maximum possible value, ki(ki — l)/2, i.e., 
Ci = 2ei/ki{ki — 1). The clustering coefficient of the whole network is the average of all 
individual C[s. 
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Figure 3. (a) An example of a primary quadrilateral, where the three outer nodes 
are directly connected to node i. (b) An example of a secondary quadrilateral, where 
one outer node (empty square) is a second neighbor of node i. 



Since there are no triangles in the hierarchial lattices, the clustering coefficient of 
every node and their average value in H(q,t) are both zero by definition. However, 
over the years generalized clustering coefficients probing higher-order loops have been 
proposed [551 EE] ■ Clearly in these hierarchial lattices the number of squares (loops of 
length 4) is significantly large, below we will seek to quantify this. 

2.5. Grid Coefficient 

As pointed out above, for hierarchial lattices the usual clustering coefficient is unable 
to quantify the order underlying their structure, which is represented by a grid-like 
frame, that can be quantified by evaluating the frequency of rectangular loops (cycles of 
length 4). We introduce the grid coefficient that allows us to uncover the presence of a 
surprising level of triangular grid ordering in the hierarchial lattices. For simplicity, we 
call cycles of length 4 quadrilaterals. The grid coefficient [56] Gi of node i is defined as 
the ratio of number of existing quadrilaterals passing by node i, Xi, to all the possible 
number of quadrilaterals attached to node i, Yj. 

Note that each quadrilateral involving node i is consist of i itself plus three outer 
nodes, according to whose nature quadrilaterals can be classified, then the grid coefficient 
can be further decomposed into two cases (see Fig. [3]): If all the outer nodes are directly 
attached to i, they form a primary quadrilateral; otherwise, if one of the outer nodes is 
a second neighbor of i, the cycle they form is a secondary quadrilateral. If a node i with 
deg : 2nd second neighbors, then the maximum possible number of primary 

quadrilaterals is Yf = 3x ( 3 i ) = ki(k i — l)(k i — 2)/2, while the maximum possible number 
of secondary quadrilaterals is Yf = ki^ n dki(ki — l)/2. In this way, for investigating the 
grid properties of the hierarchial lattices, one can define three quantities: the primary 
grid coefficient, G\ = Xf/Y?, the secondary grid coefficient G| = X?/Yf, and the total 
grid coefficient Gi = (Xf + X*)/(Y? + Yf), where Xf and Xf are the actually existing 
number of primary and secondary quadrilaterals involving i, respectively. Averaging 
these quantities over all nodes, we can obtain the respective average grid coefficients. 
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There are only secondary quadrilaterals in HLs, so G\ = 0. The analytical 
expression for the secondary grid coefficient Gl of the individual node i with degree 
k can be derived exactly. In H(q,t), for a node with degree great than two, there is a 
one-to-one correspondence between the second neighbors k 2n d of the node and its degree 
k: k2nd = -■ On the other hand, by construction, for a node with degree k > 2, for 
each of its second neighbor, there are just (J) secondary quadrilaterals passing by the 
node and the second neighbor simultaneously, thus, the existing number of secondary 
quadrilaterals is X s = -(I). Therefor, for a node of degree k > 2, the exact value of its 
grid coefficient is 

So the grid coefficient is a function of degree k, following a power-law behavior of the 
form k~ 2 for large k. It is interesting to note that a similar scaling has been observed 
in several real-life networks |56j . 



2.6. Fractality 

In fact, the hierarchial lattices grow as a inverse renormalization procedure. To find 
the fractal dimension, we follow the mathematical framework proposed in Ref . [20J . By 
construction, for large t, the different quantities grow as: 

Nt — 2q Nf-i, 

k i {t) = qk i {t-l), (19) 
Dt = 2D (W ). 

The first equation is analogous to the multiplicative process naturally found in many 
population growth systems. The second relation denotes the preferential attachment 
mechanism [6], which yields the power law degree distribution. The third equation 
describes the change of the diameter D f of the hierarchial lattices H(q,t), where D t is 
defined as by the longest shortest path between all pairs of nodes in H(q,t). 

From the relations given by Eq. ( JT9l) . we know that these quantities Nt, ki(t) and 
Df increase by a factor of 2q, q and 2, respectively. Then between any two times t±, t2 
(t\ < t 2 ), we can easily obtain the following relation: 

D i2 =2* 2 -' 1 D tlJ 

N t2 = (2qY^ N tl , (20) 
h(t 2 ) = q 1 ^ ki{h). 

From Eq. (1201) . we can derive the scaling exponents in terms of the microscopic 
parameters: the fractal dimension is ds = lr f^ , and the degree exponent of boxes 
is dk = . The exponent of the degree distribution satisfies 7 = l + ^- = 2 + |^|, 
giving the same 7 as that obtained in the direct calculation of the degree distribution, 
see Eq. QTT}. 

Note that in a class of deterministic models called pseudo-fractals, although the 
number of their nodes increase exponentially, the additive growth in the diameter with 
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time implies that these networks are small world. These models do not capture the 
fractal topology found in diverse complex networks [25l [261 12ZI [28l [291 EQl EH [32l [331 

HUES]. 

2. 1. Degree correlations 

As the field of complex networks has progressed, degree correlations [HI EJ [TQl [Til [121 
[TBI Ej have been the subject of particular interest, because they can give rise to some 
interesting network structure effects. Degree correlations can be conveniently measured 
by means of the conditional probability P(k'\k), being defined as the probability that a 
link from a node of degree k points to a node of degree k! . In uncorrelated networks, this 
conditional probability does not depend on k, it takes the form P(k'\k) = k'P(k')/(k) 



Although degree correlations are formally characterized by P(k'\k), the direct 
evaluation of the conditional probability P(k'\k) in real-life systems is a very difficult 
task, and usually gives extremely noisy results because of their finite size. To overcome 
this problem, another interesting quantity related to two node correlations, called 
average nearest-neighbor degree (ANND), has been proposed. It is a function of node 
degree, and is more convenient and practical in characterizing degree-degree correlations, 
defined by [9] 



If there are no degree correlations, Eq. ( f2~TT) gives k nn (k^ — (k^ / (fc), i.e. k nn {k^j is 
independent of k. Degree correlations are usually quantified by reporting the numerical 
value of the slope of k nn {k) as a function of node degree k. 

Degree correlations quantified by ANND have led to a first classification of complex 
networks. When k nn (k) increases with k, it means that nodes have a tendency to connect 
to nodes with a similar or larger degree. In this case the network is defined as assortative 
[TT1 [T2] . In contrast, if k nn (k) is decreasing with k, which implies that nodes of large 
degrees are likely to have the nearest neighbors with small degrees, then the network is 
said to be disassortative. 

We can exactly calculate k nn (k) for the hierarchial lattices. By construction, for 
nodes with degree greater than 2, the degrees of their neighbors are 2. Then we have 



[13j 




(21) 



k' 



k nn (k > 2) = 2. 



(22) 



For those nodes having degree 2, their average nearest-neighbor degrees are 




t<=i-l 




for q = 2 



for q > 2. 



(23) 
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Thus k nn (2) grows linearly or exponentially with time for q = 2 and q > 2, respectively. 
As the nodes with degree 2 are only connected to higher degree nodes, k nn (2) is 
significantly high. 

Degree correlations can also be described by a Pearson correlation coefficient r of 
degrees at either end of a link. It is defined as [HI [12], [33J [57] 
(fc)(Pfc nn (fc))-(fc 2 ) 2 

(k)(k 3 ) - (k 2 ) 2 ' 1 J 

If the network is uncorrelated, the correlation coefficient equals zero. Disassortative 
networks have r < 0, while assortative graphs have a value of r > 0. Substituting Eqs. 
dHJ), ( fl5l) and ( fT6l) into Eq. (Ell), we can easily see that for t > 1, r of H(q,t) is always 
negative, indicating disassortativity. 

Disassortative features in protein interaction networks were found and explained 
by Maslov and Sneppen [8] on the level of interacting proteins and genetic regulatory 
interactions. According to their results links between highly connected nodes are 
systematically suppressed, while those between highly connected and low-connected 
pairs of proteins are favored. 



2.8. Average path length 

Shortest paths play an important role both in the transport and communication within 
a network and in the characterization of the internal structure of the network. We 
represent all the shortest path lengths of H(q, t) as a matrix in which the entry dij 
is the geodesic path from node i to node j, where geodesic path is one of the paths 
connecting two nodes with minimum length. The maximum value of d^ is called the 
diameter of the network. A measure of the typical separation between two nodes in the 
hierarchical lattices is given by the average path length d t , also known as characteristic 
path length, defined as the mean of geodesic lengths over all couples of nodes at the tth 
level. 

For general q, it is not easy to derive a closed formula for the average path length 
d t of H(q,t). However, in the Appendix, we have obtained exact analytic expressions 
for d t of H(3,t), while the exact value of d t of H(2,t) has been obviously obtained in 
Ref. [50]. For q = 3 we find 



- _ 1243 ■ 2* - 475 • 2* +2 3* + 275 • 2*+ 3 3* + 275 • 2*3 2m + 19 • 2 3 * +2 3 2m - 11 • 3' +2 4 m 

* ~ 44(2 + 2'3*+ 1 )(7 + 2'3*+ 1 ) 
,iq 

r=T^2*, (25) 

leading to an exponential growth in the APL. Since in this case, N t ~ 6* for large t, 
we have d t ~ N t ogfi • While in another special case q = 2, d t ~ nI' 2 [50J. So for 
small q, the hierarchical lattices are not small worlds. We conjecture that for H(q,t), 
their APL scales as d t ~ N t ° g(2q) , which is similar to that of a hypercubic lattice 
of dimension ; — - — k. Low-dimensional regular lattices do not show the small- world 

lo g(2q) 2 ^ 
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Figure 4. (Color online) Construction of the small-world hierarchical lattices, with 
(a), (b) denoting SH 1 (2,2) and SH 2 (2,S), respectively. 

behavior of typical node-node distances. It is straightforward to show that for a regular 
lattice in D dimensions which has the shape of a square or (hyper)cube of side Z, and 
therefore has N = l D nodes, the APL increases as I, or equivalently as N l / D |58j . 

So we have shown that d t of H (q, t) has the power-law scaling behavior of the 
number of nodes N t . It is not hard to understand. As an example, let us look at the 
scheme of the growth of a particular case q = 2. Each next step in the growth of H (2, t) 
doubles the APL between a fixed pair of nodes. The total numbers of nodes and edges 
increase four-fold (asymptotically, in the infinite limit of t), see Eq. ([5]). Thus the APL 
dt of H(2, t) grows as a square power of the node number in H(2, t). 

3. Small- world hierarchical lattices 

In this section, we will discuss the construction and properties of small-world 
hierarchical lattices (SWHLs). Our goal is to reduce the diameter enough so as to 
get a logarithmically growing diameter, while maintaining the original structure of 
hierarchical lattices studied in preceding section. All these can be attained by adding a 
new central point and connecting it with a certain set of original nodes, which is akin 
to the ideas presented in Refs. [591 EH EE] ■ 

We denote the small- world hierarchical lattices as SH m (q,t). From the previous 
section, we can know that for every t > 2 and m = 1, 2, 3, • • • , t — 1, H (q, t) can be seen 
as (2qY~ m copies of H (q, m), with some node identifications. SH m (q, t) (1 < m < t — 1) 
is the graph obtained by joining a new node to every hub (node of highest degree) of 
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every copy of H(q,m), see Fig. HI In other words, the new node is connected to those 
old nodes introduced at step t — m or earlier, thus the number of new edges is exactly 
the number of vertices of H(q, t — m). For m = 0, the new graph is out of the scope of 
SH m (q, t). In this case, for simplicity, we also denote the new network SH m (q, t), where 
the central node connects to all the nodes in H(q,t). 

The diameter of this new graph depends on the value of m. We will show that, for 
some of the values of m, the diameter of SH m (q, t) exhibits a slow (logarithmic) increase 
with the total number of network nodes. Thus, this construction gives us small-world 
graphs. Next we give the properties of small-world hierarchical lattices. 

The order (number of all nodes) of SH m (q,t) is one plus the order of H(q,t). The 
size (number of all edges) of SH m (q,t) is the size of H(q,t), plus the number of added 
edges. Since the number of added edges is the order of H(q,t — m), according to Eqs. ((Sj) 
and (jSJ), we can easily see that the order and size of SH m (q } t) is N m;t = g ( 2g ) +5g-3 an( j 

= (2 g )^-(2 g )«+,(2 g )t-" + 3 g -2 > respectively . 

Because the addition of the new node has little effect on the degree distribution, 
SH m (q,t) also follow power law degree distribution with the same degree exponent 7 
as H(q,t). Additionally, for any m > 1, SH m (q,t) have no triangles, the clustering 
coefficient is zero as their counterparts H(q,t). Analogously, SH m (q,t) are self-similar 
with the identical fractal dimension d B as H(q,t) [20| |6T]. 

Different from H(q,t), SH m (q,t) have small-world property. For the sake of 
convenient expression, let us denote by Diam[H(q,t)] the diameter of H(q,t) and 
by Diam[SH m (q,t)] the diameter of SH m (q,t). Obviously, Diam[H (q, t)) = 2*. To 
compute Diam[SH m (q,t)} we need only observe that, in H(q,m), every node is at 
distance at most £j^E!iMLM!M — 2 m ~ 1 from the set of vertices of the hubs. An upper 
bound for Diam[SH m (q, t)} is Diam[H(q, m)] + 2 = 2 m + 2. It can easily be seen that 
this is also a lower bound. Therefore, Diam[SH m (q,t)} = 2 m + 2. Since the order N m j 
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of SH m (q,t) is 9( ' 2 ' ? 2 g -i' ? ~ 3 ' ^ m — 1°§2^ then Diam[SH m (q, t)] < t + 2 is small and 
scales logarithmically with the number of network nodes. Here we do not give the exact 
expression for the average path length (APL) of SH m (q,t), instead in Fig. Owe present 
APL of SH m (q,t) as a function of m. It is shown that APL becomes larger as m is 
increased. 

To summarize, in the section, we proposed a construction of small-world hierarchical 
lattices. In the construction of these small-world lattices, the underlying structure of 
the original lattices is preserved. We have shown that all of these new graphs are fractal 
and have a logarithmic diameter. 

4. Relative robustness to international attacks 

As discussed in previous section, close to many real-life networks, the hierarchical lattices 
are simultaneously self-similar and scale-free. Therefore, it is worthwhile to investigate 
the processes taking place upon them and directly compare these results with just scale- 
free networks (like BA networks). These comparisons may give us deep insight into the 
dynamic properties of networks. In the following we will investigate intentional damage 
(attack) and synchronization, respectively. This section is devoted to the robustness, 
while next section is concerned with collective synchronization behavior. 

Robustness refers to the ability of a network to avoid malfunctioning when a fraction 
of its constituents is damaged. This is a topic of obvious practical reasons, as it affects 
directly the efficiency of any process running on top of the network, and it is one of 
the first issues to be explored in the literature on complex networks [62] . Here we shall 
focus only on the topological aspects (especially self-similarity) of robustness, caused by 
targeted removal of nodes, because there is a strong correlation between robustness and 
network topology [S21 E31 ESI E51 E5! • 

One of the most important measures of the robustness of a network is its integrity, 
which is characterized by the presence of its giant connected component [62] . We call a 
network robust if it contains a giant cluster comprised of most of the nodes even after 
a fraction of its nodes are removed. Then, to know network robustness, first of all, one 
must study the variation of the giant component. 

For the study of attack vulnerability of the network, the selection procedure of the 
order in which vertices are removed is an open choice [66J . One may of course maximize 
the destructive effect at any fixed number of removed vertices. However, this requires 
the knowledge of the whole network structure, and pinpointing the vertex to attack 
in this way makes a very time-demanding computation. A more tractable choice is to 
select the vertices in the descending order of degrees in the initial network and then to 
remove vertices one by one starting from the vertex with the highest degree; this attack 
strategy will be used in the present paper. 

Figure [6] shows the performances of BA and hierarchial lattices under intentional 
attack. We plot the relative size of the largest cluster, s, after removing a fraction p 
of the largest hubs for both networks. One can find that the non-fractal scale-free BA 
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Figure 6. (Color online) Vulnerability under intentional attack of a BA network with 
average degree 4 and a H (2, 8). Both networks have the same degree exponent (7 = 3), 
the same number of nodes (43,692), and their clustering coefficient is small. Moreover, 
by construction, SH m (2, 8) have similar robustness as H(2, 8), which rule out the effect 
of average path length on targeted attack. Thus, the difference in the resilience seen 
in this figure is attributed to fractality and the different degree of anticorrelations. 

networks are more sensitive to sabotage of a small fraction of the nodes, leading support 
to the view of Song et. al. [20]. While both networks collapse at a finite fraction p c , 
evidenced by the decrease of s toward zero, the fractal network has a significantly larger 
threshold (p c ~ 0.004) compared to the non-fractal threshold (p c m 0.001), the former 
threshold is about 4 times than the latter, suggesting a significantly higher robustness 
of the fractal networks to intentional attacks. Also, it is interesting to note that for the 
hierarchical lattices, s is a function of p with a staircase-like form. 

It is not strange at all that a giant connected component in self-similar networks is 
robust against the targeted deletion of nodes, while non-fractal scale-free networks are 
extremely vulnerable to targeted attacks on the hub. In non-fractal topologies, the hubs 
are connected and form a central compact core, such that the removal of a few of the 
largest hubs has catastrophic consequences for the network. In self-similar networks, 
hubs are more dispersed (see Fig. |5J), their disassortativity and self-similar property 
significantly increases the robustness against targeted attacks. This could explain why 
some real-life networks have evolved into a fractal and disassortative architecture [19, 20J. 

5. Synchronization 

The ultimate goal of the study of network structure is to study and understand the 
workings of systems built upon those networks. Recently, along with the study of 
purely structural and evolutionary properties [TJ [2], there has been increasing interest 
in the interplay between the dynamics and the structure of complex networks 0, H] . One 
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particular issue attracting much attention is the synchronizability of oscillator coupling 
networks [67]. Synchronization is observed in diverse natural and man-made systems 
and is directly related to many specific problems in a variety of different disciplines. 
It has found practical applications in many fields including communications, optics, 
neural networks and geophysics [HEl EH [TO, ED [721 E3]- After studying the relevant 
characteristics of network structure, which is described in the previous sections, we will 
study the synchronization behavior on the networks. 

We follow the general framework proposed in [7^1 175] , where a criterion based on 
spectral techniques was established to determine the stability of synchronized states 
on networks. Consider a network of N identical dynamical systems with linearly and 
symmetric coupling between oscillators. The set of equations of motion for the system 
are 



where Xj = F(xj) governs the dynamics of each individual node, H(xj) is the output 
function and a the coupling strength, and Gij is the Laplacian matrix, defined by 
Gu = ki if the degree of node i is k iy = — 1 if nodes % and j are connected, and 
G^ = otherwise. 

Since matrix G is positive semidefinite and each rows of it has zero sum, all 
eigenvalues of G are real and non-negative and the smallest one is always equal to zero. 
We order the eigenvalues as = Ai < A2 < • • ■ < Xn- Then one can use the ratio of the 
maximum eigenvalue A at to the smallest nonzero one A2 to measure the synchronizability 
of the network [7_U [75]. If the eigenratio R = X^/X 2 satisfies R < a 2 /ai, we say the 
network is synchronizable. Here the eigenratio R depends on the the network topology, 
while a 2 /«i depends exclusively on the dynamics of individual oscillator and the output 
function. Ratio R = Xn/X 2 represents the synchronizability of the network: the larger 
the ratio, the more difficult it is to synchronize the oscillators, and vice versa. 

After reducing the issue of synchronizability to finding eigenvalues of the Laplacian 
matrix G, we now investigate the synchronization of our networks. Figure [7J shows the 
synchronizabily of H(2,t), SH m (2,5), SH m (2,6), as well as the BA networks. One can 
see that for the same network order, R of the BA networks is much smaller than that 
of hierarchical lattices H(2,t), which implies that the synchronizability of the former is 
much better. While for SH m (2, 5) and SH m (2, 6), the dependence relation of eigenratio 
R(m) on m is more complicated: R{0) < R{2) < R(l) < R(x\x > 2); for m > 2, R(m) 
increases with m > 2. 

Why coupling systems on the BA networks, H(2,t) and SH m (2,t) exhibit very 
different synchronizability? Previously reported results have indicated that underlying 
network structures play significant roles in the synchronizability of coupled oscillators. 
However, the key structural feature that determines the collective synchronization 
behavior remains unclear. Many works have discussed this issue. Some authors believe 
that shorter APL tends to enhance synchronization [711 ESI [77] . In contrast, Nishikawa 
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Figure 7. (Color online) (Top) The eigenratio R as a function of network order TV for 
BA networks with average degree 4 and H(2,t). All quantities for BA networks are 
averaged over 50 realizations. (Bottom) The dependence of eigenratio R of SH m (2,5) 
and SH m (2, 6) on the value of m. 



et al. reported that synchronizability is suppressed as the degree distribution becomes 
more heterogeneous, even for shorter APL [78]. In Ref. [79], the authors asserted that 
larger average node degree corresponds to better synchroizability. 

All these may rationally explain the relations between synchroizability and network 
structure in some cases, but do not well account for the difference of synchroizability 
between the graphs under consideration: as known from the preceding section, for fixed 
t, the APL of SH m (2,t) increases with m, but R does not always decrease with m; 
BA networks and H (2, t) have identical degree exponent 7, while their synchroizability 
differs very much; in addition, the average degree of SH\(2,t) is higher than that 
of SH 2 (2,t), but the former is more difficult to synchronize than the latter. All 
these show that the degree distribution is generally not sufficient to characterize 
the synchronizability of scale-free networks [80, 81], and that smaller average path 
length does not necessarily predict better synchroizability [81]. We speculate that the 
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synchroizability on BA networks is better than on H(2,t) and SH m (2,t) rests mainly 
with the self-similar structure. The genuine reasons need further research. 

6. Conclusions 

In conclusion, we have studied a family of deterministic networks called hierarchical 
lattices (HLs) from the perspective of complex networks. The deterministic self-similar 
construction allow one to derive analytic exact expressions for the relevant features of 
HLs. Our results shows that HLs exhibit many interesting topological properties: they 
follow a power-law degree distribution with exponent tuned from 2 to 3; their clustering 
coefficient is null but their grid coefficient follows a power-law phenomenon; they are not 
small-world, the APL scales like a power-law in the number of nodes; they have a fractal 
topology with a general fractal dimension; and they are disassortative networks. Our 
results indicate it is not true that a power-law degree alone create small-world networks 
[82] . and further support the conjecture that scale-free networks with fractal scaling 
are disassortative. The disassortativity property is ease to understand by checking the 
growth process of HLs, where the rich (large nodes) get richer but at the expense of the 
poor (small nodes). In other words, the hubs prefer to grow by connections with less- 
connected nodes rather to other hubs, which leads to disassortativity. So we have found 
a good example — hierarchical lattices — which show that self-similar scale-free networks 
are preferably disassortative in their degree-degree correlations. 

We have also introduced a deterministic construction of a family of small-world 
hierarchical lattices (SWHLs) and investigated their topological characteristics. We 
have shown that the some basic structure features of the hierarchical lattices (HLs) are 
preserved, including degree distribution, clustering coefficient and fractality, while the 
small-world phenomenon arises. 

In addition, we have studied the dynamical processes such as intention damage 
and collective synchronization and have shown the comparisons between HLs and BA 
networks as well as SWHLs. We have found that self-similarity and disassortativity 
increase the robustness of networks under intentional attacks. Some qualitative 
explanations have been given showing that a fractal and disassortative topology 
structure is more robust. Although HLs and SWHLs have relative better robustness, 
they exhibit poorer synchronizability than BA networks based possibly on the same 
reasons as that of their insensitiveness to sabotage. We have shown that sychronizability 
of scale-free networks is not an intrinsic property of the exponent of degree distribution, 
and that small APL does not imply good synchronizability. 

In future, it would be worth searching for other stochastic networks displaying 
finite fractal dimension spectra. Moreover, it is more interesting that the presence of 
self-similarity and disassortativity, as well as the absence of small-world properties of 
the HLs might have a relevant effect on the other dynamic process such as epidemic 
spreading [83], routing traffic [EUE5], games [86] and transport [87] taking place on the 
networks. 
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Appendix A. Derivation of the average path length for q = 3 

We denote the set of nodes constituting the hierarchical lattices H(q,t) after t 
construction steps as L q j. Then the APL for L qjt is defined as: 

d t = , — -, , (A.l) 

N t (N t - I) /2 1 ; 

where 

D t = Yl d v ( A ' 2 ) 

denotes the sum of the chemical distances between two nodes over all pairs, and d^ is 
the chemical distance between nodes i and j. Although there are some difficulties in 
obtaining a closed formula for d t holding true for all q, the hierarchical lattices have a 
self-similar structure that allows one to calculate d t analytically according to different 
q. As shown in Fig. \A1\ the lattice L q t+1 may be obtained by joining 2q copies of L q t at 

-(a) 
J q,t ! 

as 



the hubs, which are labeled as 1$ , a = 1, 2, • • • , 2q. Then we can write the sum D t +i 



D t+1 = 2qD t + A t} (A.3) 

where A t is the sum over all shortest paths whose endpoints are not in the same L qjt 
branch. The solution of Eq. IA.3I is 

t-i 

D t = (2q) t - 1 D 1 + ^(2g)'- m - 1 A m . (A.4) 

771=1 

The paths that contribute to A t must all go through at least one of the q + 2 edge 
nodes (for q = 3 see Fig. IA11 where A, B, C, D, E are the five edge nodes) at which 
the different L q t branches are connected. The analytical expression for A t for general 
q, called the crossing paths, is not easy to derive. We trace the formula only for the 
particular case of q = 3 as follow. 

In what follows, we write L^ t as L t for brevity. Denote A"" 8 as the sum of all 
shortest paths with endpoints in and L t . If and uf meet at an edge node, 



A"' p rules out the paths where either endpoint is that shared edge node. If L\ and 

do not meet, A) 
the total sum A t is 



do not meet, A"" 5 excludes the paths where either endpoint is any edge node. Then 



A, = A]' 2 + A, 1 ' 3 + A, 1 ' 4 + A, 1 ' 5 + A]' 6 + A 2 / 
+ A 2 / + A 2 ' 5 + A 2 ' 6 + A 3 / + A 3 ' 5 + A 3 / 

+ A 4 ' 5 + A 4 ' 6 + A J 5 ' 6 - 5 • 2 t+1 . (A.5) 
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The last term at the end compensates for the overcounting of certain paths: the shortest 
path between A and B, with length 2' +1 , is included in A*' 6 , A 2 ' 5 and A 3 ' 4 ; the shortest 
path between C and E, with length 2* +1 , is included in both A*' 3 and A 4 ' 6 ; the shortest 
path between D and E, with length 2 i+1 , is included in both A 2 ' 3 and A 4 ' 5 ; the shortest 
path between C and D, also with length 2* +1 , is included in both A*' 2 and A^' 6 . 

By symmetry, A, 1 ' 2 = A]' 3 = A 2 ' 3 = A t 5 ' 6 = A 4,5 = A 4,6 = A t 1,6 = A 2,5 = A 3,4 and 

a 1,4 A 1,5 a 2,4 a 2.6 a 3,5 a 3,6 ,i . 

A t ' — A t ' — A t ' = A t ' — A t ' — A t ' , so that 

A« = 9A]' 2 + 6Aj' 4 - 5 • 2 m , (A.6) 
where A*' 2 is given by the sum 



A t 2 = Yl di i 

i€L^\j&Lf\i^A,j^A 

22 (d iA + d A j) 

- z L?\i^A,&A 

= (N t - 1) d iA + (N t -i)Yl d ^ 

= 2(JV t -l) d iA , (A.7) 

where ^2 ieL (r) d iA = X^ 6 l< 2 ) ^4? have been used. To find 5^ ieL (i) diA, we examine the 

structure of the hierarchical lattice at the tth level. In \ there is v t {m) points with 
diA — m i where 1 < m < 2', and z/ 4 (m) can be written recursively as 

{3* if m is odd 

,m . (A.8) 

v t-i\ — ) 11 m 1S even 

We can write y}. c7 -(i) in terms of vArn) as 

2* 

Eqs. (IA.8|) and fl A. 9|) relate f t and /t_i, which allow one to resolve f t by induction as 
follow: 

2 t-i 2 t-i 

/ t = ^)(2A;-l)3 t + ^;2A;.i/ t _ 1 (A;) 
fc=i fc=i 
= 3*2 2 *~ 2 + 2/t_i 

= I 2 *- 2 (l4 + 6 m ), (A.10) 
5 

where / x = ^i(l) + 2^(2) = 5 has been used. Substituting Eq. flATTUl) and iV t = z± f^ 
into Eq. f]A.7[) . we obtain 

A]' 2 = — 2*- 1 (3 x 6* + 2)(6 m + 14) . (A.ll) 
25 
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E 




C 

Figure Al. For q — 3, the hierarchical lattice after t + 1 construction steps, L t +i, is 
composed of six copies of L t denoted as (\ — 1, 2, • • • , 6), which are connected to 
one another as above. 

Continue analogously, 

ieL™ , jeL™ ,&A,C, j^B,E 

J2 + 2 * + d iz) 

i£L { t\ j£L ( t 4) ,i^A, j^E, d iA +d jE <2 t 

+ Yl ( d * c + 2 ' + d ^ 

i&L^, jeL ( t 4 \i^C, jjtB, d tC +d jB <2t 

+ J2 2 * +1 - ( A - 12 ) 

i&L^, je4 A \ i^A, j+E, d lA +d jE =2* 

The first terms equal the second ones and are denoted by g t , and the third term is 
denoted by h t , so that Aj' = 2g t + h t . One can compute the quantity g t as 

2*-2 2*-l-m 

9t = ^ v t {m)v t {m!){m + 2* + mf) 

m=l m'=l 

2'- 1 -2 2'- 1 -l-fc 

vt-i(k)vt-i(k')(2k + 2 t + 2k l ) 

k=l k'=l 
2 t-l_i 2 t - 1 -k 

+ vt-i(k)3 t (2k + 2 t + 2k' - 1) 

k=l k'=l 
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2 t-l_x 2 t - 1 -k 



Yl 3 t u t _ 1 (k')(2k- 1 + 2* + 2k') 



(A.13) 



fc=i k'=i 

+ ^ 3 2 *(2A; - 1 + 2* + 2fc' - 1) . 

k=l k'=l 

The fourth terms can be summed directly, yielding 

3 2 *2*- 3 (2* - 2) 2 + 3 2 * _1 2*- 1 (2* _1 + 1)(2* - 2). (A.14) 

In Eq. ( 1A.13I) . the second and third terms are equal to each other and can be simplified 
by first summing over k', yielding 

2^-1 

3* £ ^(A;)^ 2 '- 2 -^-*; 2 ). 



(A.15) 



k=l 



For use in Eq. ([QKjl . ^Li ft-i(*0 = #t-i - 2, and using Eq. (TA~T0l . 

^ A;^_ 1 (A;) = ^A;^ 1 (A:)-2*- 1 



fc=i 



fc=i 



(A.16) 



(A.17) 



(A.18) 



= 2*- d (6*-6)/5. 
Similarly to Eq. flA.lOj) . we get 

y fcVi(fc) = l^-^*- 1 - -e*- 1 + — 4*- 1 . 

^ 1 u ; 5 2 10 

fc=i 

With these results, Eq. ( 1A.15I) becomes 

^3*(2 • 4*- 1 6 < - 1 - 3 • 4*- 1 + 6*" 1 ) . 

With Eqs. flATil) and (jATT8|) . Eq. flA~T3l becomes 

1 „+ 3 
? 

(A.19) 

Considering the initial condition g\ = 0, we can solve Eq. ( 1A.19I) inductively leading to 

_ A 2 *-3(_i7i + 17 . 2 *+ 2 3* +1 + 5 ■ 2* +3 3 2 * +1 
y 85 V 

-85-9*-17-4* +1 9*). (A.20) 



9t = 2g t _ 1 + -3*6*4* - -3*6* - -3*4* + 2*~ 3 4*9* 
6 6 ' 

- -4*9* + ^2*9* . 
2 2 



To find an expression for A 1 / : , now the only thing left is to evaluate h t as 



2*-l 



tH = 2* +1 ^z/ t (m)^(2*-m) 

m=l 

= 2^|> 2 (m) 



m=l 
2 t-l 



)t+l 



5>« + £ iSt-i(k) 



k=l 



k=l 



16 2 * + 2/l t _! 



(A.21) 
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where we have used the the symmetry v t [m) = v t {2 1 — m). Since h\ = 36, Eq. (1A.21D is 
solved inductively: 

^ = ^2' +1 (18'-l). (A.22) 

From Eqs. flA~20l) and flA~22l . 

A. 1 - 4 = - A 2 *- 2 (-m + 17 • 2*+ 2 3 m + 5 • 2* +3 3 2 ' +1 - 85 • 9* - 17 • 4 m 9*) 
* 85 V ; 

+ ^2 m (18'-l). (A.23) 

Substituting Eqs. (1A.11I) and (1A.23I) into Eq. (1A.6D . we obtain the final expression for 
the crossing paths A t : 

A t = —2 t - 1 (2 + 3 ■ 6*)(14 + 6 m ) - — 2*- 1 (-171 + 17 • 2* +2 3 m + 5 ■ 2* +3 3 2m 
25 85 

54 

- 85 ■ 9* - 17 • 4 m 9*) + ^2 m (18* - 1) - 5 ■ 2 t+1 , (A.24) 
Substituting Eqs. (1A.24I) for A m into Eq. (jA.4[) . and using D\ = 14, we have 

D t = — !— (1243 ■ 2* - 475 ■ 2 t+2 3* + 275 • 2* +3 3* + 275 ■ 2*3 2t+1 
2200 v 

+ 19 • 2 3t+2 3 2t+1 - 11 • 3 t+2 4 m ) (A.25) 

Inserting Eq. (IA.25I) into Eq. (jA.ip . one can obtain the analytical expression for d t in 
Eq. ([25]). 
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Abstract. In this paper, firstly, we study analytically the topological features of a 
family of hierarchical lattices (HLs) from the view point of complex networks. We 
derive some basic properties of HLs controlled by a parameter q: scale-free degree 
distribution with exponent 7 = 2 + null clustering coefficient, power-law behavior 
of grid coefficient, exponential growth of average path length (non-small- world), fractal 
scaling with dimension ds = ^wf^ , and disassortativity. Our results show that 
scale-free networks are not always small-world, and support the conjecture that self- 
similar scale-free networks are not assortative. Secondly, we define a deterministic 
family of graphs called small- world hierarchical lattices (SWHLs). Our construction 
preserves the structure of hierarchical lattices, including its degree distribution, fractal 
architecture, clustering coefficient, while the small-world phenomenon arises. Finally, 
the dynamical processes of intentional attacks and collective synchronization arc 
studied and the comparisons between HLs and Barabasi- Albert (BA) networks as well 
as SWHLs are shown. We find that the self-similar property of HLs and SWHLs 
significantly increases the robustness of such networks against targeted damage on 
hubs, as compared to the very vulnerable non fractal BA networks, and that HLs have 
poorer synchronizability than their counterparts SWHLs and BA networks. We show 
that degree distribution of scale-free networks does not suffice to characterize their 
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synchronizability, and that networks with smaller average path length are not always 
easier to synchronize. 



PACS numbers: 89.75.Da, 05.45.Df, 36.40.Qv, 05.45.Xt 
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1. Introduction 

Topological characteristics, such as scale-free degree distribution, small-world effect, 
fractal scaling and degree correlations, have recently attracted much attention in 
network science. The last few years have witnessed a tremendous activity devoted 
to the characterization and understanding of networked systems [U EJ [3j Hj. Small- 
world property [5] and scale-free behavior [6j are two unifying concepts constituting 
our basic understanding of the organization of real-life complex systems. Small-world 
property refers to the one that the expected number of edges (links) needed to pass 
from one arbitrarily selected node (vertex) to another one is low, which grows at most 
logarithmically with the number of nodes. Scale-free behavior means the majority of 
nodes in a network have only a few connections to other nodes, whereas some nodes are 
connected to many other nodes in the network. This poses a fundamental question how 
these two characteristics are related. It has been observed that small-world property 
and scale-free behavior are not independent p]: scale-free networks, normally, have 
extremely short average path length (APL), scaling logarithmically or slower with 
system size. Is this universal? 

In fact, the above mentioned two properties (i.e. small-world property and scale- 
free behavior) do not provide sufficient characterizations of the real-world systems. It 
has been observed that real networks exhibit ubiquitous degree correlations among 
their nodes [H [9j [TOj EH 121 E31 ITS] . This translates in the observation that the 
degrees of nearest neighbor nodes are not statistically independent but mutually 
correlated in practically every network imaginable. Correlations play an important 
role in the characterization of network topology, and have led to a first classification 
of complex networks A series of recent measurements indicate social networks 

are all assortative, while all biological and technological networks disassortative: in 
social networks there is a tendency for the hubs to be linked together, in biological and 
technological network the hubs show the opposite tendency, being primarily connected 
to less connected nodes. Correlations are now a very relevant issue, especially in view 
of the important consequences that they can have on dynamical processes taking place 
on networks [151 ESI EH EE] • 

Recently, it has been discovered by the application of a renormalization procedure 
that diverse real networks, such as the WWW, protein-protein interaction networks 
and metabolic networks, exhibit fractal scaling and topological self-similarity [T9J, EH1 
I2TI [22| [23] . Fractal scaling implies that in a network the minimum number of node- 
covering boxes Nb of linear size £b scales with respect to ^ as a power-law Nb, with an 
exponent that is given by a finite fractal dimension ds [19] . Self-similarity refers to the 
invariant scale-free distribution probability to find a node with degree k, P(k) ~ fc -7 , 
i.e. the exponent 7 remains the same under the renormalization with different box 
sizes [T9l [23] . In complex networks, fractality and self-similarity do not always imply 
each other: A fractal network model is self-similar, while a self-similar network is not 
always fractal [23]. One can obtain the fractal dimension ds by measuring the ratio of 
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Nb over the total number of nodes N in the network, which satisfies Nb/N ~ Iff 11 1 . After 
renormalizing the networks, the degree /cb(^b) of a node in the renormalized network 
versus the largest degree khub inside the box that was contracted to one node with degree 
khub in the renormalization process exhibits a scaling behavior: kg (&b) = s(£B)kh u b, 
where s(£b) is assumed to scales like s(£b) ~ £~B dk with dk being the degree exponent of 
the boxes. In self-similar scale-free networks, the relation between the three indexes 7, 
ds and dk satisfies 7 = 1 + d B /d k [19|. 

Correlations and topological fractality are important properties for many real- 
world complex systems. Then a natural fundamental question is raised how the 
two characteristics relate to each other. Recent researches [201 122] have shown that 
self-similar scale-free networks are not assortative, and the qualitative feature of 
disassortativity is scale-invariant under renormalization. Moreover, self-similarity and 
disassortativity of scale-free networks make such networks more robust against a sinister 
attack on nodes with large degree, as compared to the very vulnerable non fractal scale- 
free networks [2D] . 

However, do small-world property and scale-free behavior always go along? How 
do systems have evolved into self-similar disassortative scale-free networks? How the 
dynamical processes such as intentional attack and synchronization are influenced by 
the topological fractality and disassortativity of scale-free networks? Such a series of 
important questions still remain open. To relate these questions, in this paper we 
therefore launch a study seeking a better understanding of the relations among these 
topological properties. 

It is of interest to study above important questions with deterministic methods. 
Because of their strong advantages, deterministic network models have received much 
attention [251 12H1 [2TI [28l [291 EQl EU [32l [331 EH [35l EHl EZl EH [39] . First, the method of 
generating deterministic networks makes it easier to gain a visual understanding of how 
networks are shaped, and how do different nodes relate to each other [25]; moreover, 
deterministic networks allow to compute analytically their topological properties, which 
have played a significant role, both in terms of explicit results and a guide to and a 
test of simulated and approximate methods I251I2S1I2ZII2E1I2S1ED1EI1E21E31E31ES1 
I3P1 [37] [38] [3D] . On the other hand, deterministic networks can be easily extended to 
produce random variants which exhibit the classical characteristics of many real-life 
systems 0D1 SB S2J M, IB] • 

Inspired by the above mentioned questions, here we first introduce a deterministic 
family of networks. These networks are called hierarchical lattices (HLs), which yield 
exact renormalization-group solutions [45j |46l HTJ EHJ |49l [50]. From the perspective 
of complex network, we show that HLs are simultaneously scale-free, self-similar and 
disassortative, but lack the small-world property. Then we present a deterministic 
construction of a class of small- world hierarchical lattices (SWHLs), which preserve 
the basic structure properties including power law degree distribution, self-similarity, 
and disassortativeness, while lead to the small- world effect. Finally, we investigate the 
effects of network structures on the dynamics taking place in them. 
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Figure 1. Iterative construction method of the hierarchical lattices for some limiting 
cases. 

This article is organized as follows: In Sec. II, we introduce the construction of 
hierarchical lattices (HLs) and study their topological features including the degree 
distribution, moments, clustering coefficient, grid coefficient, self-similarity, degree 
correlations, and average path length (APL). The detailed exact derivation about APL 
is shown in Appendix. In Sec. Ill, we propose the deterministic construction of the 
small-world hierarchical lattices (SWHLs) and study their properties. In Sec. IV, attack 
tolerance of HLs is studied and the comparisons between HLs and Barabasi- Albert (BA) 
networks are shown. In Sec. V, we do a comparative investigation of synchronization in 
HLs, SWHLs and BA networks. Sec. VI is devoted to our conclusions. 

2. Hierarchical lattices 

In this section, from topological perspective of complex networks, we present the 
construction and the basic properties such as degree distribution, clustering coefficient, 
average path length (APL), fractality, and correlations of the hierarchical lattices (HLs). 

2.1. Construction of the Lattice 

The hierarchical lattices [H] are constructed in an iterative manner as shown in Fig. [TJ 
We denote the hierarchical lattices (networks) after t generations by H(q,t), q>2 and 
t > 0. The networks are constructed as follows: For t — 0, H (q, 0) is an edge connecting 
two points. For t > 1, H(q,t) is obtained from H(q,t — 1). We replace each of the 
existing edges in H(q,t — 1) by the connected cluster of edges on the right of Fig. [U 
The growing process is repeated t times, with the infinite lattices obtained in the limit 
t — > oo. Figure [2] shows the growing process of the networks for three particular cases 
of q = 2, q = 3, and q — 4. It should be noted that in the hierarchical lattice of q = 2 
case [15] , the Migdal-Kadanoff J5TJ [52] recursion relations with dimension 2 and length 
rescaling factor 2 are exact. 

Griffiths and Kaufman provided two explanations for the construction of the 
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a) 

q=2 




t=0 



t=1 



t=2 



Figure 2. Examples of the hierarchical lattices for some particular case of q = 2, 
q = 3, and q = 4, showing the first three steps of the iterative process. 



hierarchical lattices [37], which are called "aggregation" and "miniaturization". In 
essence, these two interpretations reflect the self-similar structure of the hierarchical 
lattices, which allow one to calculate analytically their topological characteristics. 

Next we compute the numbers of nodes (vertices) and links (edges) in H(q,t). Let 
L v {t) and L e (t) be the numbers of vertices and edges created at step t, respectively. 
Note that each of the existing edges yields q nodes, and the addition of each new node 
leads to two new edges. By construction, for t > 1, we have 

L v {t) = qL e {t - 1) (1) 

and 



L e (t) = 2L v (t). 

Considering the initial condition L v (0) = 2 and L e (0) = 1, it follows that 
L v (t)=q(2qy- 1 

and 

L e (t) = {2qf. 

Thus the number of total nodes N t and edges E t present at step t is 



(2) 
(3) 
(4) 

(5) 
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and 

Et = L e (t) = (2qf, (6) 

respectively. 

2.2. Degree distribution 

Let ki(t) be the degree of node i at step t. Then by construction, it is not difficult to 
find following relation: 

ki(t) = qki(t-l), (7) 

which expresses a preference attachment [Sj. If node i is added to the network at step 
ti, ki(ti) = 2 and hence 

hit) = 2^. (8) 

Therefore, the degree spectrum of the network is discrete. It follows that the degree 
distribution is given by 

L v (0) 2 



P(k) 



Nt - g(2g)«+3g-2 f ° r *i - 

L V (U)_ qjqh fQTt . >1 (9) 

2q-l 

otherwise 



and that the cumulative degree distribution [HI [26] is 

p (k) - V Lv{p) - q{2q)U + 39 " 2 (10) 

In * 

Substituting for ^ in this expression using ti — t— gives 

_ g(2g) i (|)"^+3g-2 
cuml j ~ g(2g)* + 3q - 2 

for large t. (11) 



, _f1 + ln2 



2 

So the degree distribution follows a behavior power law with the exponent 7 = 2 + 



ln2 
In g ' 

For q = 2, Eq. ( TTTT) recovers the result of the particular case p = previously obtained 
in Ref. [26] . 



£5. Moments 

Information on how the degree is distributed among the nodes of a undirected network 
can be obtained either by the degree distribution P(k), or by the calculation of the 
moments of the distribution. The n-moment of P(k) is defined as: 

(k n ) = J2k n P(k). (12) 
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The first moment (k) is the mean degree. At arbitrary step t, the average vertex 
degree of H(q, t) is 

2E t 2(2g-l)(2 g )' 

<*>' = 74 = g (2,)- + 3,-2 ' (13) 

For large t, it is small and approximately equal to a finite value 4 — -. 

We can also calculate higher moments of the distribution P(k). For instance, the 
second moment, which measures the fluctuations of the connectivity distribution, is 
given by 

(k 2 )t = ^Y,n v {U)[k{t u t)}\ (14) 

* u=o 

where k(ti, t) is the degree of a node at step t which was generated at step £j. This 
quality expresses the average of degree square over all nodes in the network. It has large 
impact on the behavior of dynamical processes taking place in networks 
Substituting Eqs. d2J), ((SD and (jSJ) into Eq. ffT4|) . we derive 

4g(2g - 1) q* - 2* 2(2g - l)g* 



^ 4 f (3t + 3) 
4* + 2 
4(g-l)(2g-l) ,q 



for q > 2 
for g = 2 

— . . i for g > 2 . , 

g(g-2) \2) H (15) 

3(t + l) forg = 2 

In this way, second moment of degree distribution (k 2 ) has been calculated explicitly, 
and the result shows that it becomes infinite for large t. In fact, because the degree 
exponent 7 < 3, all n-moments (n > 2) diverge. For example, we can analogously get 
the third moment as 

2*+2^3t+3 _ 2 t + 1 q 3t + 2 _|_ 2*+ 3 g 3 *+ 1 _ 2 t+2 q 3t — 2 2t+A q t+1 + 2 2t+3 q t 
\ k )* = 2 2 *g i+3 - 2 2t+1 q t+1 + 3 ■ 2*g 3 - 3g ■ 2* +1 - 2* +1 g 2 + 2 t + 2 

For the special case of g = 2, it reduces to 
q . 9 3t _ a . ? 2t 

= 22t+ \ , (17) 
which diverges as an exponential law when t is very large. 



2.4- Clustering Coefficient 

The clustering coefficient defines a measure of the level of cohesiveness around any given 
node. By definition, the clustering coefficient |5] Cj of node i is the ratio between the 
number of edges e« that actually exist among the ki neighbors of node i (i.e. the number 
of triangles attached to a vertex %) and its maximum possible value, ki(ki — l)/2, i.e., 
Ci = 2ei/ki{ki — 1). The clustering coefficient of the whole network is the average of all 
individual C[s. 



Self- similarity, small-world, scale-free scaling, disassortativity, and robustness in hierarchical lattices^ 



a) W I 



* - - ■+ 




Figure 3. (a) An example of a primary quadrilateral, where the three outer nodes 
are directly connected to node i. (b) An example of a secondary quadrilateral, where 
one outer node (empty square) is a second neighbor of node i. 



Since there are no triangles in the hierarchial lattices, the clustering coefficient of 
every node and their average value in H(q,t) are both zero by definition. However, 
over the years generalized clustering coefficients probing higher-order loops have been 
proposed [551 EE] ■ Clearly in these hierarchial lattices the number of squares (loops of 
length 4) is significantly large, below we will seek to quantify this. 

2.5. Grid Coefficient 

As pointed out above, for hierarchial lattices the usual clustering coefficient is unable 
to quantify the order underlying their structure, which is represented by a grid-like 
frame, that can be quantified by evaluating the frequency of rectangular loops (cycles of 
length 4). We introduce the grid coefficient that allows us to uncover the presence of a 
surprising level of triangular grid ordering in the hierarchial lattices. For simplicity, we 
call cycles of length 4 quadrilaterals. The grid coefficient [56] Gi of node i is defined as 
the ratio of number of existing quadrilaterals passing by node i, Xi, to all the possible 
number of quadrilaterals attached to node i, Yj. 

Note that each quadrilateral involving node i is consist of i itself plus three outer 
nodes, according to whose nature quadrilaterals can be classified, then the grid coefficient 
can be further decomposed into two cases (see Fig. [3]): If all the outer nodes are directly 
attached to i, they form a primary quadrilateral; otherwise, if one of the outer nodes is 
a second neighbor of i, the cycle they form is a secondary quadrilateral. If a node i with 
deg : 2nd second neighbors, then the maximum possible number of primary 

quadrilaterals is Yf = 3x ( 3 i ) = ki(k i — l)(k i — 2)/2, while the maximum possible number 
of secondary quadrilaterals is Yf = ki^ n dki(ki — l)/2. In this way, for investigating the 
grid properties of the hierarchial lattices, one can define three quantities: the primary 
grid coefficient, G\ = Xf/Y?, the secondary grid coefficient G| = X?/Yf, and the total 
grid coefficient Gi = (Xf + X*)/(Y? + Yf), where Xf and Xf are the actually existing 
number of primary and secondary quadrilaterals involving i, respectively. Averaging 
these quantities over all nodes, we can obtain the respective average grid coefficients. 
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There are only secondary quadrilaterals in HLs, so G\ = 0. The analytical 
expression for the secondary grid coefficient Gl of the individual node i with degree 
k can be derived exactly. In H(q,t), for a node with degree great than two, there is a 
one-to-one correspondence between the second neighbors k 2n d of the node and its degree 
k: k2nd = -■ On the other hand, by construction, for a node with degree k > 2, for 
each of its second neighbor, there are just (J) secondary quadrilaterals passing by the 
node and the second neighbor simultaneously, thus, the existing number of secondary 
quadrilaterals is X s = -(I). Therefor, for a node of degree k > 2, the exact value of its 
grid coefficient is 

So the grid coefficient is a function of degree k, following a power-law behavior of the 
form k~ 2 for large k. It is interesting to note that a similar scaling has been observed 
in several real-life networks |56j . 



2.6. Fractality 

In fact, the hierarchial lattices grow as a inverse renormalization procedure. To find 
the fractal dimension, we follow the mathematical framework proposed in Ref . [20J . By 
construction, for large t, the different quantities grow as: 

Nt — 2q Nf-i, 

k i {t) = qk i {t-l), (19) 
Dt = 2D (W ). 

The first equation is analogous to the multiplicative process naturally found in many 
population growth systems. The second relation denotes the preferential attachment 
mechanism [6], which yields the power law degree distribution. The third equation 
describes the change of the diameter D f of the hierarchial lattices H(q,t), where D t is 
defined as by the longest shortest path between all pairs of nodes in H(q,t). 

From the relations given by Eq. ( JT9l) . we know that these quantities Nt, ki(t) and 
Df increase by a factor of 2q, q and 2, respectively. Then between any two times t±, t2 
(t\ < t 2 ), we can easily obtain the following relation: 

D i2 =2* 2 -' 1 D tlJ 

N t2 = (2qY^ N tl , (20) 
h(t 2 ) = q 1 ^ ki{h). 

From Eq. (1201) . we can derive the scaling exponents in terms of the microscopic 
parameters: the fractal dimension is ds = lr f^ , and the degree exponent of boxes 
is dk = . The exponent of the degree distribution satisfies 7 = l + ^- = 2 + |^|, 
giving the same 7 as that obtained in the direct calculation of the degree distribution, 
see Eq. QTT}. 

Note that in a class of deterministic models called pseudo-fractals, although the 
number of their nodes increase exponentially, the additive growth in the diameter with 
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time implies that these networks are small world. These models do not capture the 
fractal topology found in diverse complex networks [25l [261 12ZI [28l [291 EQl EH [32l [331 

HUES]. 

2. 1. Degree correlations 

As the field of complex networks has progressed, degree correlations [HI EJ [TQl [Til [121 
[TBI Ej have been the subject of particular interest, because they can give rise to some 
interesting network structure effects. Degree correlations can be conveniently measured 
by means of the conditional probability P(k'\k), being defined as the probability that a 
link from a node of degree k points to a node of degree k! . In uncorrelated networks, this 
conditional probability does not depend on k, it takes the form P(k'\k) = k'P(k')/(k) 



Although degree correlations are formally characterized by P(k'\k), the direct 
evaluation of the conditional probability P(k'\k) in real-life systems is a very difficult 
task, and usually gives extremely noisy results because of their finite size. To overcome 
this problem, another interesting quantity related to two node correlations, called 
average nearest-neighbor degree (ANND), has been proposed. It is a function of node 
degree, and is more convenient and practical in characterizing degree-degree correlations, 
defined by [9] 



If there are no degree correlations, Eq. ( f2~TT) gives k nn (k^ — (k^ / (fc), i.e. k nn {k^j is 
independent of k. Degree correlations are usually quantified by reporting the numerical 
value of the slope of k nn {k) as a function of node degree k. 

Degree correlations quantified by ANND have led to a first classification of complex 
networks. When k nn (k) increases with k, it means that nodes have a tendency to connect 
to nodes with a similar or larger degree. In this case the network is defined as assortative 
[TT1 [T2] . In contrast, if k nn (k) is decreasing with k, which implies that nodes of large 
degrees are likely to have the nearest neighbors with small degrees, then the network is 
said to be disassortative. 

We can exactly calculate k nn (k) for the hierarchial lattices. By construction, for 
nodes with degree greater than 2, the degrees of their neighbors are 2. Then we have 



[13j 




(21) 



k' 



k nn (k > 2) = 2. 



(22) 



For those nodes having degree 2, their average nearest-neighbor degrees are 




t<=i-l 




for q = 2 



for q > 2. 



(23) 
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Thus k nn (2) grows linearly or exponentially with time for q = 2 and q > 2, respectively. 
As the nodes with degree 2 are only connected to higher degree nodes, k nn (2) is 
significantly high. 

Degree correlations can also be described by a Pearson correlation coefficient r of 
degrees at either end of a link. It is defined as [HI [12], [33J [57] 
(fc)(Pfc nn (fc))-(fc 2 ) 2 

(k)(k 3 ) - (k 2 ) 2 ' 1 J 

If the network is uncorrelated, the correlation coefficient equals zero. Disassortative 
networks have r < 0, while assortative graphs have a value of r > 0. Substituting Eqs. 
dHJ), ( fl5l) and ( fT6l) into Eq. (Ell), we can easily see that for t > 1, r of H(q,t) is always 
negative, indicating disassortativity. 

Disassortative features in protein interaction networks were found and explained 
by Maslov and Sneppen [8] on the level of interacting proteins and genetic regulatory 
interactions. According to their results links between highly connected nodes are 
systematically suppressed, while those between highly connected and low-connected 
pairs of proteins are favored. 



2.8. Average path length 

Shortest paths play an important role both in the transport and communication within 
a network and in the characterization of the internal structure of the network. We 
represent all the shortest path lengths of H(q, t) as a matrix in which the entry dij 
is the geodesic path from node i to node j, where geodesic path is one of the paths 
connecting two nodes with minimum length. The maximum value of d^ is called the 
diameter of the network. A measure of the typical separation between two nodes in the 
hierarchical lattices is given by the average path length d t , also known as characteristic 
path length, defined as the mean of geodesic lengths over all couples of nodes at the tth 
level. 

For general q, it is not easy to derive a closed formula for the average path length 
d t of H(q,t). However, in the Appendix, we have obtained exact analytic expressions 
for d t of H(3,t), while the exact value of d t of H(2,t) has been obviously obtained in 
Ref. [50]. For q = 3 we find 



- _ 1243 ■ 2* - 475 • 2* +2 3* + 275 • 2*+ 3 3* + 275 • 2*3 2m + 19 • 2 3 * +2 3 2m - 11 • 3' +2 4 m 

* ~ 44(2 + 2'3*+ 1 )(7 + 2'3*+ 1 ) 
,iq 

r=T^2*, (25) 

leading to an exponential growth in the APL. Since in this case, N t ~ 6* for large t, 
we have d t ~ N t ogfi • While in another special case q = 2, d t ~ nI' 2 [50J. So for 
small q, the hierarchical lattices are not small worlds. We conjecture that for H(q,t), 
their APL scales as d t ~ N t ° g(2q) , which is similar to that of a hypercubic lattice 
of dimension ; — - — k. Low-dimensional regular lattices do not show the small- world 

lo g(2q) 2 ^ 
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behavior of typical node-node distances. It is straightforward to show that for a regular 
lattice in D dimensions which has the shape of a square or (hyper) cube of side Z, and 
therefore has N = l D nodes, the APL increases as I, or equivalently as iV 1//D [58] . 

So we have shown that d t of H(q,t) has the power-law scaling behavior of the 
number of nodes N t . It is not hard to understand. As an example, let us look at the 
scheme of the growth of a particular case q = 2. Each next step in the growth of H (2, t) 
doubles the APL between a fixed pair of nodes. The total numbers of nodes and edges 
increase four- fold (asymptotically, in the infinite limit of t), see Eq. ([5]). Thus the APL 
d t of H (2, t) grows as a square power of the node number in H(2, t). 

3. Small- world hierarchical lattices 

In this section, we will discuss the construction and properties of small-world 
hierarchical lattices (SWHLs). Our goal is to reduce the diameter enough so as to 
get a logarithmically growing diameter, while maintaining the original structure of 
hierarchical lattices studied in preceding section. All these can be attained by adding a 
new central point and connecting it with a certain set of original nodes, which is akin 
to the ideas presented in Refs. [59j EH El] . 

We denote the small- world hierarchical lattices as SH m (q,t). As mentioned in the 
previous section, for every t > 2 and m = 1, 2, 3, • • • , t — 1, H(q, t) can be seen as (2g)' _m 
copies of H(q,m), with some node identifications. SH m (q,t) (1 < m < t — 1) is the 
graph obtained by joining a new node to every hub (node of highest degree) of every 
copy of H(q,m), see Fig. [2j In other words, the new node is connected to those old 
nodes introduced at step t — m or earlier, thus the number of new edges is exactly the 
number of vertices of H(q,t — m). For m = 0, the new graph is out of the scope of 
SH m (q, t). In this case, for simplicity, we also denote the new network SH m (q, t), where 
the central node connects to all the nodes in H(q,t). 

The diameter of this new graph depends on the value of m. We will show that, for 
some of the values of m, the diameter of SH m (q, t) exhibits a slow (logarithmic) increase 
with the total number of network nodes. Thus, this construction gives us small-world 
graphs. Next we give the properties of small-world hierarchical lattices. 

The order (number of all nodes) of SH m (q,t) is one plus the order of H(q,t). The 
size (number of all edges) of SH m (q,t) is the size of H(q,t), plus the number of added 
edges. Since the number of added edges is the order of H(q, t — m), according to Eqs. [5] 
and El we can easily see that the order and size of SH m (q,t) is N m j = q ( 2q ^+^ q ~ 3 anc l 

E m>t = (2g)t+1 - (2g)t 2 ;! ( f )t ' m+3g - 2 , respectively. 

Because the addition of the new node has little effect on the degree distribution, 
SH m (q,t) also follow power law degree distribution with the same degree exponent 7 
as H(q,t). Additionally, for any m > 1, SH m (q,t) have no triangles, the clustering 
coefficient is zero as their counterparts H(q,t). Analogously, SH m (q,t) are self-similar 
with the identical fractal dimension d B as H(q,t) [201 EI]- We denote the small-world 
hierarchical lattices as SH m (q,t). As mentioned in the previous section, for every t > 2 
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Figure 4. (Color online) Construction of the small-world hierarchical lattices, with 
(a), (b) denoting SH 1 (2,2) and SH 2 (2,S), respectively. 

and m = 1, 2, 3, • • • , t — 1, H(q, t) can be seen as {2q) t ~ m copies of H (q, m), with some 
node identifications. SH m (q,t) (1 < m < t — 1) is the graph obtained by joining a new 
node to every hub (node of highest degree) of every copy of H(q,m), see Fig. HI In 
other words, the new node is connected to those old nodes introduced at step t — m or 
earlier, thus the number of new edges is exactly the number of vertices of H(q,t — m). 
For m = 0, the new graph is out of the scope of SH m (q,t). In this case, for simplicity, 
we also denote the new network SH m (q,t), where the central node connects to all the 
nodes in H(q, t). 

The diameter of this new graph depends on the value of m. We will show that, for 
some of the values of m, the diameter of SH m (q, t) exhibits a slow (logarithmic) increase 
with the total number of network nodes. Thus, this construction gives us small-world 
graphs. Next we give the properties of small-world hierarchical lattices. 

The order (number of all nodes) of SH m (q,t) is one plus the order of H(q,t). The 
size (number of all edges) of SH m (q,t) is the size of H(q,t), plus the number of added 
edges. Since the number of added edges is the order of H(q, t — m), according to Eqs. [5] 
and El we can easily see that the order and size of SH m (q,t) is iV m>i = g ( 2g ^ ^ g ~ 3 and 
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E m>t = (2^-(2 g )^-" + 3 g -2 > respectiyely . 

Because the addition of the new node has little effect on the degree distribution, 
SH m (q,t) also follow power law degree distribution with the same degree exponent 7 
as H(q,t). Additionally, for any m > 1, SH m (q,t) have no triangles, the clustering 
coefficient is zero as their counterparts H(q,t). Analogously, SH m (q,t) are self-similar 
with the identical fractal dimension ds as H(q,t) [20| |6T]. 

Different from H(q,t), SH m (q,t) have small- world property. For the sake of 
convenient expression, let us denote by Diam[H (q, t)] the diameter of H(q,t) and 
by Diam[SH m (q,t)] the diameter of SH m (q,t). Obviously, Diam[H(q,t)] = 2*. To 
compute Diam[SH m (q,t)] we need only observe that, in H(q,m), every node is at 
distance at most Dmm Wi< m )] — 2 m ~ 1 from the set of vertices of the hubs. An upper 
bound for Diam[SH m (q,t)] is Diam[H(q,m)} + 2 = 2 m + 2. It can easily be seen that 
this is also a lower bound. Therefore, Diam[SH m (q, t)] = 2 m + 2. Since the order N m> t of 
SH m (q, t) is g(M-±j£zj? ; if 772 < log 2 1, then Diam[SH m (q, t)] < t + 2 is small and scales 
logarithmically with the number of network nodes. We do not give the exact expression 
for the average path length (APL) of SH m (q,t), instead in Fig. [5] we present APL of 
SH m (q,t) as a function of m. It is shown that APL becomes larger as m is increased. 

To summarize, in the section, we proposed a construction of small-world hierarchical 
lattices. In the construction of these small-world lattices, the underlying structure of 
the original lattices is preserved. We have shown that all of these new graphs are fractal 
and have a logarithmic diameter. 

4. Relative robustness to international attacks 

As discussed in previous section, close to many real-life networks, hierarchical-lattices 
are simultaneously self-similar and scale-free. Therefore, it is worthwhile to investigate 
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the processes taking place upon hierarchical-lattices and directly compare these results 
with just scale-free networks (like BA networks). These comparisons may give us deep 
insight into the dynamic properties of networks. In the following we will investigate 
intentional damage (attack) and synchronization, respectively. This section is devoted to 
the robustness, while next section is concerned with collective synchronization behavior. 

Robustness refers to the ability of a network to avoid malfunctioning when a fraction 
of its constituents is damaged. This is a topic of obvious practical reasons, as it affects 
directly the efficiency of any process running on top of the network, and it is one of 
the first issues to be explored in the literature on complex networks [62]. Here we shall 
focus only on the topological aspects (especially self-similarity) of robustness, caused by 
targeted removal of nodes, because there is a strong correlation between robustness and 
network topology [E21 M, Hfl M, M\ ■ 

One of the most important measures of the robustness of a network is its integrity, 
which is characterized by the presence of its giant connected component [62] . We call a 
network robust if it contains a giant cluster comprised of most of the nodes even after 
a fraction of its nodes are removed. Then, to know network robustness, first of all, one 
must study the variation of the giant component. 

For the study of attack vulnerability of the network, the selection procedure of the 
order in which vertices are removed is an open choice [66J . One may of course maximize 
the destructive effect at any fixed number of removed vertices. However, this requires 
the knowledge of the whole network structure, and pinpointing the vertex to attack 
in this way makes a very time-demanding computation. A more tractable choice is to 
select the vertices in the descending order of degrees in the initial network and then to 
remove vertices one by one starting from the vertex with the highest degree; this attack 
strategy will be used in the present paper. 

Figure [6] shows the performances of BA and hierarchial lattices under intentional 
attack. We plot the relative size of the largest cluster, s, after removing a fraction p 
of the largest hubs for both networks. One can find that the non-fractal scale-free BA 
networks are more sensitive to sabotage of a small fraction of the nodes, leading support 
to the view of Song et. al. [20]. While both networks collapse at a finite fraction p c , 
evidenced by the decrease of s toward zero, the fractal network has a significantly larger 
threshold (p c ~ 0.004) compared to the non-fractal threshold (p c ~ 0.001), the former 
threshold is about 4 times than the latter, suggesting a significantly higher robustness 
of the fractal networks to intentional attacks. Also, it is interesting to note that for the 
hierarchical lattices, s is a function of p with a staircase-like form. 

It is not strange at all that a giant connected component in self-similar networks is 
robust against the targeted deletion of nodes, while non-fractal scale-free networks are 
extremely vulnerable to targeted attacks on the hub. In non-fractal topologies, the hubs 
are connected and form a central compact core, such that the removal of a few of the 
largest hubs has catastrophic consequences for the network. In self-similar networks, 
hubs are more dispersed (see Fig. [2]), their disassortativity and self-similar property 
significantly increases the robustness against targeted attacks. This could explain why 
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Figure 6. (Color online) Vulnerability under intentional attack of BA network with 
average degree 4 and a H(2, 8). Both networks have the same degree exponent (7 = 3), 
the same number of nodes (43,692), and their clustering coefficient is small. Moreover, 
by construction, SH m (2, 8) have similar robustness as H(2, 8), which rule out the effect 
of average path length on targeted attack. Thus, the difference in the resilience seen 
in this figure is attributed to fractality and the different degree of anticorrelations. 



some real-life networks have evolved into a fractal and disassortative architecture [19, 20J. 
5. Synchronization 

The ultimate goal of the study of network structure is to study and understand the 
workings of systems built upon those networks. Recently, along with the study of 
purely structural and evolutionary properties [U [2], there has been increasing interest 
in the interplay between the dynamics and the structure of complex networks [3], 0] . One 
particular issue attracting much attention is the synchronizability of oscillator coupling 
networks [67]. Synchronization is observed in diverse natural and man-made systems 
and is directly related to many specific problems in a variety of different disciplines. 
It has found practical applications in many fields including communications, optics, 
neural networks and geophysics [Ml EH [TQl [TTJ [721 After studying the relevant 
characteristics of network structure, which is described in the previous sections, we will 
study the synchronization behavior on the networks. 

We follow the general framework proposed in [74l [75] , where a criterion based on 
spectral techniques was established to determine the stability of synchronized states 
on networks. Consider a network of iV identical dynamical systems with linearly and 
symmetric coupling between oscillators. The set of equations of motion for the system 
are 



where Xj = F(x,) governs the dynamics of each individual node, H(xj) is the output 



N 




(26) 
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function and o the coupling strength, and G^ is the Laplacian matrix, defined by 
Gu = ki if the degree of node i is fcj, Gij = — 1 if nodes % and j are connected, and 
G^ = otherwise. 

Since matrix G is positive semidefinite and each rows of it has zero sum, all 
eigenvalues of G are real and non-negative and the smallest one is always equal to zero. 
We order the eigenvalues as = Ai < A 2 < • ■ ■ < Xn- Then one can use the ratio of the 
maximum eigenvalue A at to the smallest nonzero one A 2 to measure the synchronizability 
of the network [7H [75]. If the eigenratio R = \n/^2 satisfies R < a^/ai, we say the 
network is synchronizable. Here the eigenratio R depends on the the network topology, 
while 0:2/0:1 depends exclusively on the dynamics of individual oscillator and the output 
function. Ratio R = represents the synchronizability of the network: the larger 

the ratio, the more difficult it is to synchronize the oscillators, and vice versa. 

After reducing the issue of synchronizability to finding eigenvalues of the Laplacian 
matrix G, we now investigate the synchronization of our networks. Figure [7J shows the 
synchronizabily of H(2,t), SH m (2,5), SH m (2,6), as well as the BA networks. One can 
see that for the same network order, R of the BA networks is much smaller than that 
of hierarchical lattices H(2,t), which implies that the synchronizability of the former is 
much better. While for SH m {2, 5) and SH m {2, 6), the dependence relations of eigenratio 
R(m) on m more complicated: -R(O) < R{2) < R(l) < R(x\x > 2); for m > 2, R(m) 
increases with m > 2. 

Why coupling systems on the BA networks, H(2,t) and SH m (2,t) exhibit very 
different synchronizability? Previously reported results have indicated that underlying 
network structures play significant roles in the synchronizability of coupled oscillators. 
However, the key structural feature that determines the collective synchronization 
behavior remains unclear. Many works have discussed this issue. Some authors believe 
that shorter APL tends to enhance synchronization [7U [76J [77J . In contrast, Nishikawa 
et al. reported that synchronizability is suppressed as the degree distribution becomes 
more heterogeneous, even for shorter APL [78]. In Ref. [79], the authors asserted that 
larger average node degree corresponds to better synchroizability. 

All these may rationally explain the relations between synchroizability and network 
structure in some cases, but do not well account for the difference of synchroizability 
between the graphs under consideration: as known from the preceding section, for fixed 
t, the APL of SH m (2,t) increases with m, but R does not always decrease with m; 
BA networks and H (2, t) have identical degree exponent 7, while their synchroizability 
differs very much; in addition, the average degree of SHi(2,t) is higher than that 
of S H2(2, t), but the former is more difficult to synchronize than the latter. All 
these show that the degree distribution is generally not sufficient to characterize 
the synchronizability of scale-free networks [80j [81], and that smaller average path 
length does not necessarily predict better synchroizability [81] . We speculate that the 
synchroizability on BA networks is better than on H(2,t) and SH m (2,t) rests mainly 
with the self-similar structure. The genuine reasons need further research. 
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Figure 7. (Color online) (Top) The eigenratio R as a function of network order TV for 
BA networks with average degree 4 and H(2,t). All quantities for BA networks are 
averaged over 50 realizations. (Bottom) The dependence of eigenratio R of SH m (q, 5) 
and SH m (q, 6) on the value of m. 



6. Conclusions 

In conclusion, we have studied a family of deterministic networks called hierarchical 
lattices (HLs) from the perspective of complex networks. The deterministic self-similar 
construction allow one to derive analytic exact expressions for the relevant features of 
HLs. Our results shows that HLs exhibit many interesting topological properties: they 
follow a power-law degree distribution with exponent tuned from 2 to 3; their clustering 
coefficient is null but their grid coefficient follows a power-law phenomenon; they are not 
small-world, the APL scales like a power-law in the number of nodes; they have a fractal 
topology with a general fractal dimension; and they are disassortative networks. Our 
results indicate it is not true that a power-law degree alone create small-world networks 
[82J, and further support the conjecture that scale-free networks with fractal scaling 
are disassortative. The disassortativity property is ease to understand by checking the 
growth process of HLs, where the rich (large nodes) get richer but at the expense of the 
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poor (small nodes). In other words, the hubs prefer to grow by connections with less- 
connected nodes rather to other hubs, which leads to disassortativity. So we have found 
a good example — hierarchical lattices — which show that self-similar scale-free networks 
are preferably disassortative in their degree-degree correlations. 

We have also introduced a deterministic construction of a family of small-world 
hierarchical lattices (SWHLs) and investigated their topological characteristics. We 
have shown that the some basic structure features of the hierarchical lattices (HLs) are 
preserved, including degree distribution, clustering coefficient and fractality, while the 
small-world phenomenon arises. 

In addition, we have studied the dynamical processes such as intention damage 
and collective synchronization and have shown the comparisons between HLs and BA 
networks as well as SWHLs. We have found that self-similarity and disassortativity 
increase the robustness of networks under intentional attacks. Some qualitative 
explanations have been given showing that a fractal and disassortative topology 
structure is more robust. Although HLs and SWHLs have relative better robustness, 
they exhibit poorer synchronizability than BA networks based possibly on the same 
reasons as that of their insensitiveness to sabotage. We have shown that sychronizability 
of scale-free networks is not an intrinsic property of the exponent of degree distribution, 
and that small APL does not imply good synchronizability. 

In future, it would be worth searching for other stochastic networks displaying 
finite fractal dimension spectra. Moreover, it is more interesting that the presence of 
self-similarity and disassortativity, as well as the absence of small-world properties of 
the HLs might have a relevant effect on the other dynamic process such as epidemic 
spreading [83], routing traffic [HI [85], games [86] and transport [87] taking place on the 
networks. 
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Appendix A. Derivation of the average path length for q = 3 

We denote the set of nodes constituting the hierarchical lattices H(q,t) after t 
construction steps as L Qjt - Then the APL for L Q)t is defined as: 



D t 



(A.l) 



N t (N t - l)/2 



where 



D, 



't = 





(A.2) 



denotes the sum of the chemical distances between two nodes over all pairs, and dij is 
the chemical distance between nodes % and j. Although there are some difficulties in 
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obtaining a closed formula for d t holding true for all q, the hierarchical lattices have a 
self-similar structure that allows one to calculate d t analytically according to different 
q. As shown in Fig. IA11 the lattice L qtt +i may be obtained by joining 2q copies of L q>t at 
the hubs, which are labeled 1, 2, • • • , 2q. Then we can write the sum D t+1 

as 

D t+1 = 2qD t + A t , (A.3) 

where A t is the sum over all shortest paths whose endpoints are not in the same L qjt 
branch. The solution of Eq. IA.3I is 

t-i 

D t = (2q) t - 1 D 1 + ^(2g)'- m - 1 A m . (A.4) 

771=1 

The paths that contribute to A t must all go through at least one of the q + 2 edge 
nodes (for q = 3 see Fig. IA11 where A, B, C, D, E are the five edge nodes) at which 
the different L q ^ branches are connected. The analytical expression for At for general 
q, called the crossing paths, is not easy to derive. We trace the formula only for the 
particular case of q = 3 as follow. 

In what follows, we write L 3 t as L t for brevity. Denote Aj as the sum of all 
shortest paths with endpoints in l{ and L t . If Li and Lj. meet at an edge node, 
At' 13 rules out the paths where either endpoint is that shared edge node. If L[ and 

do not meet, A"" 3 excludes the paths where either endpoint is any edge node. Then 
the total sum A t is 

A t = A 1 / + A, 1 ' 3 + A, 1 ' 4 + A, 1 ' 5 + A]' 6 + Af' 3 
+ Af + Af + A 2 ' 6 + Af + A 3 / + Af 

+ A 4 ' 5 + A 4 ' 6 + A J 5 ' 6 - 5 • 2 t+1 . (A.5) 

The last term at the end compensates for the overcounting of certain paths: the shortest 
path between A and B, with length 2' +1 , is included in Af , A 2 ' 5 and A 3 ' 4 ; the shortest 
path between C and E, with length 2* +1 , is included in both A 1 ' 3 and A 4,6 ; the shortest 
path between D and E, with length 2* +1 , is included in both Af and A 4,5 ; the shortest 
path between C and D, also with length 2' +1 , is included in both Af and Af . 

By symmetry, Af = Af = A 2,3 = Af = A 4,5 = A 4,6 = Af = Af = Af and 

a 1,4 A 1,5 A 2,4 A 2,6 a 3,5 A 3,6 ,i , 

A t = A t = A t = A t = A t = A t , so that 

A t = 9Af + 6Af - 5 ■ 2 m , (A.6) 
where Af is given by the sum 

Af = d v 

ieL?>, j€L[ 2) ,i^A,j^A 
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E 




C 

Figure Al. For q — 3, the hierarchical lattice after t + 1 construction steps, L t +i, is 
composed of six copies of L t denoted as (\ — 1, 2, • • • , 6), which are connected to 
one another as above. 



tel.™ jeL™ 
2(N t -l)J2 d iA, (A.7) 



where Y)-^tW diA = JZ-^ri 2 ) &Ai have been used. To find Y^rCO diA, we examine the 

structure of the hierarchical lattice at the tth level. In Lf\ there is z/ t (m) points with 
diA = m, where 1 < m < 2', and z/ t (m) can be written recursively as 

{3* if m is odd 

( m . (A.8) 

v t-\\—) !i m is even 

We can write J2-crW ^a in terms of v t {rn) as 

2* 

ft = rf ij4 = m • i/ t (m) . (A. 9) 

Eqs. ( ]A.8j) and flA.9j) relate and /t_i, which allow one to resolve f t by induction as 
follow: 

2 t-i 2 t-i 
/ t = ^(2fc-l)3* + J2 2k-u t ^(k) 

k=l k=l 
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= 3*2 2t ~ 2 + 2/ t _i 

= -2'- 2 (14 + 6 m ), (A.10) 
5 

where f x = + 2v x {2) = 5 has been used. Substituting Eq. (jATIOl) and N t = I± ^- 
into Eq. (jA.70 . we obtain 

A i )2 = J_ 2 *-i(3 x 6 * + 2 )(6 m + 14) . (A.ll) 
25 

Continue analogously, 

+ 2 * + ^) 

+ ( d * c + 2 * + ^ 

+ 2 t+1 . (A.12) 

The first terms equal the second ones and are denoted by gt, and the third term is 
denoted by h t , so that Aj' = 2g t + h t . One can compute the quantity g t as 

2*-2 2*-l-m 

#t = ^ v t (m)v t (m')(m + 2* + m') 

m=l m'=l 
2 *-i_2 2'- 1 -l-fc 



]T ^ ^_i(A:)^_ 1 (A;')(2fc + 2* + 2A:') 

k=i k'=i 

+ ^-i(A:)3 t (2fc + 2* + 2Jfe'- 1) 

fc=i fe'=i 

2 t-i_i 2*- 1 -fe 

+ ^t-i{k')(2k- 1 + 2* + 2fc') 



fe=i fc'=i 

2 t-l_l 2'- 1 -fc 

+ 5^ 3 2 *(2£; - 1 + 2* + 2Jfc' - 1) . (A.13) 

k=l k'=l 

The fourth terms can be summed directly, yielding 

3 2 *2*~ 3 (2* - 2) 2 + 3 2t ' 1 2 t ~ 1 (2 t - 1 + 1)(2* - 2). (A.14) 

In Eq. (1A.13j) . the second and third terms are equal to each other and can be simplified 
by first summing over k', yielding 

2^-1 

3* z/ t _i(/c)(3 • 2 2 *~ 2 - 2 l k - k 2 ) . (A. 15) 

fc=i 
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For use in Eq. flAlol) . EfcTi" 1 "t-i(*0 = ^-i ~ 2 > and usin g E <1- dSHOJ), 

^ A;^ 1 (A;) = ^^ 1 (A;)-2*- 1 
fc=i 

= 2*- 3 (6'-6)/5. (A.16) 



k=l 



Similarly to Eq. flA.lOj) . we get 



^ * u 1 5 2 10 

fe=i 



With these results, Eq. (1A.15I) becomes 
1 



3*(2 • 4*- 1 6*~ 1 - 3 • A 1 - 1 + 6* -1> 
With Eqs. flA~T4l) and flA~l8|) . Eq. flA~T3l becomes 

g t = 2g t ^ + ^3*6*4* - -3*6* - -3*4* + 2*" 3 4*9* 
6 6 4 

1 

~ 2 " " ' 2" 

Considering the initial condition = 0, we can solve Eq. (1A.19j) inductively leading to 

3 



(A.17) 



(A.18) 



-4*9* + -2*9* . 



(A.19) 



gt 



— 2*~ 3 (-171 + 17 • 2* +2 3* +1 + 5 • 2* +3 3 2 * +1 
85 

85 ■ 9* - 17- 4* +1 9*) . 



(A.20) 



To find an expression for A^' 4 , now the only thing left is to evaluate h t as 

2'-l 



/> t = 2 t+1 J> t (m> t (2* 

m=l 



m) 



m=l 

2 t-l 2 t_1 -l 



)t+l 



16 2 * + 2/^-1 



fc=l 



fc=l 



(A.21) 



where we have used the the symmetry vt(m) = z/ t (2* — m). Since fox = 36, Eq. (IA.21I) is 
solved inductively: 

9 



^ = ^2* +1 (18*-l). 
From Eqs. fTA~20l and flA~22l . 



(A.22) 



1,4 



t+3n2t+l 



2*~ 2 (-171 + 17 • 2* +2 3* +1 + 5 • 2* +3 3 

85 V 

+ A 2 *+ 1 (18 t - 1) . 



35 • 9* - 17 • 4* +i 9*) 
(A.23) 
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Substituting Eqs. OA-llO and 0A.230 into Eq. (IA.60 . we obtain the final expression for 
the crossing paths A t : 

A t = —2 t -\2 + 3 • 6*)(14 + 6 m ) - — 2*- 1 (-171 + 17 • 2* +2 3 m + 5 ■ 2* +3 3 2m 
25 85 

54 

- 85 ■ 9* - 17 • 4 m 9*) + — 2 m (18* - 1) - 5 ■ 2 t+1 , (A.24) 
Substituting Eqs. (I A. 241) for A m into Eq. (LA. 41) . and using Di = 14, we have 

D* = — (1243 • 2* - 475 • 2* +2 3* + 275 • 2* +3 3* + 275 • 2*3 2t+1 
2200 v 

+ 19 ■ 2 34+2 3 2m - 11 ■ 3 t+2 4 m ) (A.25) 

Inserting Eq. (1A.25I) into Eq. (1A.1I) . one can obtain the analytical expression for d t in 
Eq. (El. 
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